function [lagmat]=ExtractLagVec(a,lag)
lagmat=nan(length(a),lag);
for i=1:lag
    lagmat(lag+1:end,i)=a(i:length(a)-lag-1+i);
end